Non-diffracting Optical Beams in a Three-level Raman System 
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Diffractionless propagation of optical beams through atomic vapors is investigated. The atoms in 
the vapor are operated in a three-level Raman configuration. A suitably chosen control beam couples 
to one of the transitions, and thereby creates a spatially varying index of refraction modulation in the 
warm atomic vapor for a probe beam which couples to the other transition in the atoms. We show 
that a Laguerre-Gaussian control beam allows to propagate single Gaussian probe field modes as well 
as multi-Gaussian modes and non-Gaussian modes over macroscopic distances without diffraction. 
This opens perspectives for the propagation of arbitrary images through warm atomic vapors. 
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I. INTRODUCTION 

A long sought-after goal is the all-optical processing 
of information, replacing state-of-the-art electronic han- 
dling. It is well known, however, that major obstacles 
arise in particular at spatial scales small compared to 
the involved wavelengths. Light cannot be focused to 
arbitrarily small areas, and it cannot directly be used 
to image or write structures below a scale given by the 
involved wavelength [l|, Q . A key origin for these restric- 
tions is diffraction, leading, e.g., to the Rayleigh limit 
for the case of imaging [3]. Due to the enormous sig- 
nificance for applications, e.g., in nano- or life-sciences, 
a number of schemes have been invented to circumvent 
or even surpass the standard limits in a various setups. 
Examples are sub- wavelength microscopy [3-fl5j , and op- 
tical sub- wavelength lithography [16|-|22| . But more fun- 
damentally, diffraction also already impacts the propa- 
gation of a light beam in free space. For example, a 
Gaussian probe beam with a certain width at the focal 
plane spreads as it propagates away from the focal plane, 
on a length scale given by the Rayleigh length [2]. This 
paraxial diffraction spreading occurs in free space as well 
as inside a medium, and affects applications such as the 
transmission of small images, the transfer of small mask 
structures to the target surface in lithography, or the 
steering or optical manipulation of light beams. 

It was recognized that this spreading due to parax- 
ial diffraction can be decreased, increased or completely 
eliminated by propagating the light beam through a 
medium with spatially varying optical properties. This 
can be achieved via the use of external coherent control 
fields, inducing a spatially varying susceptibility and in- 
ducing phase shifts via different physical processes [23| . 
The spatial profile of a strong control beam renders the 
medium inhomogeneous along the transverse direction. 
Since absorption and dispersion of the medium are cru- 



cially dependent on the intensity and the shape of the 
strong control beam, the control beam will strongly af- 
fect the propagation of the probe beam. The ability to 
control the spatially varying susceptibility of the medium 
via external fields holds the key to eventually being able 
to control the spreading of optical beams. These and re- 
lated techniques are of significance also because they have 
the potential of complementing optical information trans- 
mission with all-optical information processing, thereby 
replacing the currently ubiquitous electronic information 
processing. 

Thus it is not surprising that a great number of 
schemes both in theory and in experiment demonstrated 
the ability to optically control the propagation of light. 
Most of these works can be understood in analogy to a 
waveguide, confining the light in specific modes transver- 
sally. This could enable one to implement app lications 
originally proposed for physical waveguides [3lL l32j | with 
ail-optically induced waveguides. For example, it was 
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FIG. 1: (Color online) Diffractionless beam propagation 
through a vapor cell, (a) shows a schematic setup to study 
the propagation of a probe laser field through an atomic vapor 
cell. The propagation dynamics is controlled using an addi- 
tional control laser field with a Laguerre-Gaussian transversal 
profile, (b) Three-level A configuration of the vapor atoms. 
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shown that spatial variations in a medium induced by a 
control laser field can lead to focusing or de-focusing of a 
probe beam [331-36] . An optically written waveguide real- 
ized by applying a Laguerre-Gaussian control field near- 
resonantly to an atomic medium to propagate a probe 
beam was demonstrated in (2J| and discussed in [25|, [26[ . 
Waveguiding in a warm rubidium vapor based on an off- 
resonant Raman configuration was reported in [27| . A 
signal beam was shown to propagate with small spot size 
over several diffraction lengths. In a variant of such se- 
tups, the spatial profile of a control beam can be cast 
onto a probe beam. Recently, this concept was imple- 
mented in Rubidium atomic vapor, based on coherent 
population trapping (28|. In this experiment, a spatial 
profile was transferred from the coupling beam onto the 
probe beam, which was found to exhibit feature sized five 
times smaller than the diffraction-limited size. 

However, all of these schemes suffer from the fact that 
the modes other than the specific supported mode can- 
not propagate without diffraction, or even not at all. 
This prohibits the desired transmission of arbitrary im- 
ages. In order to overcome this limitation, very recently, 
a scheme for the direct propagation of arbitrary probe 
beam images without diffraction was put forward theo- 
retically [29[ and verified experimentally [30| . This setup 
is based on an atomic three-level system reminiscent of 
electromagnetically induced transparency, but assumes 
an explicit breaking of the usual two-photon resonance 
condition such that strong absorption occurs. The un- 
derlying physics in this setup is fundamentally different 
from previous setups, as it relies on the movement of 
the atoms in order to compensate the diffraction, rather 
than on a stationary waveguide. A method to propagate 
optical images with reduced diffraction based on the sat- 
urated absorption mechanism in a two-level atom was 
proposed in [37|. Future ideas along this line therefore 
can be expected to advance key applications in optical in- 
formation transmission and processing, as well in reading 
and writing at the sub-wavelength scale. 

In this paper, we discuss the light propagation through 
waveguides written by a spatially structured control laser 
field into a warm atomic vapor. The Laguerre-Gaussian 
control field and the probe laser field couple with the 
atoms to form a three-level system in A configuration 
as shown in Fig. [1] We start by deriving analytical re- 
sults for the steady state susceptibility experienced by 
the probe beam. We then consider a thermal medium by 
averaging the susceptibility over a Maxwellian velocity 
distribution of the vapor atoms. The probe beam prop- 
agation dynamics is analyzed by numerically solving the 
Maxwell-Schrdinger equation including diffraction and 
dispersion in two spatial dimensions transversal to the 
propagation direction, the spatial profile of the control 
field, and the thermal motion of the vapor atoms. We 
find that for suitable parameter choices, a single Gaus- 
sian probe field mode can be propagated diffraction-less 
or even with sub-diffraction through the warm vapor. We 
then proceed to analyze more complex mode structures 



of the probe beam, and show that a Gaussian double 
pulse can be propagated without distortion or diffrac- 
tion through the optically written waveguide. Finally, we 
demonstrate that a beam formed out of multiple secant- 
hyperbolic shaped pulses can be propagated diffraction- 
less through the warm atomic vapor. The results on non- 
Gaussian beam propagation open perspectives for the 
propagation of arbitrary images in warm atomic vapors. 

The paper is organized as follows. In Sec. [Til we intro- 
duce our model, discuss the susceptibility of a thermal va- 
por, and describe the Maxwell-Schrodinger equation used 
to analyze the propagation dynamics of the probe beam 
through the optically written waveguide. In Sec. IIII Al 
we discuss the static properties of the medium suscepti- 
bility under the action of the spatially structured control 
laser beam. In Sec. IIII Bt we then discuss our results on 
the propagation of different probe beams through the op- 
tically written waveguide. Finally, Sec. IIVI discusses and 
summarizes our results. 



II. THEORETICAL CONSIDERATIONS 

A. The model 

We consider the geometry as shown in Fig. [TJ A weak 
probe laser field with frequency ui p is focused at the en- 
try of an atomic vapor, and co-propagating with an ad- 
ditional control laser field with frequency u> c focused at 
the back of the vapor cell. The cell is filled with a gas 
of three level atoms in A-configuration as shown in Fig. 
QJb) . The electric fields of the two beams can be defined 
as 

Ej{r,t) = Ef\r,t) + c.c. , (la) 
E^ + \r,t) = e,-£ 3 -(r,t) e ~ iu i t+at *' r , (lb) 

where Ej is the slowly varying envelope and e,j is unit 
polarization vector of the field. The index j £ {c,p} 
denotes the probe (p) or control (c) field, respectively. 
The polarization vectors of the beams are chosen in such 
a way that the control beam E c is tuned to the atomic 
transition |1) ■<-» |3) while the probe beam E p is coupled 
to the atomic transition |1) «-» |2). The Hamiltonian 
of the A-system under electric-dipole and rotating-wave 
approximation can then be written as 

H =H Q + Hi , (2a) 
ff =fiwi3|l><l|+fej23|2)<2|, (2b) 

H I = -{\l){2\d l2 -E^ 

+ |l)(3|d 13 -^ +) + H.c), (2c) 

where uijs denotes the resonance frequency on the tran- 
sition \j) «-» |3), d\j = (l\d\j) is the matrix element of 
the electric dipole moment operator d, and we have set 
the energy of state |3) to zero. We now make use of the 
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unitary transformation 



W = e 
U = fiw c |l)(l| + h(uj c 

to rewrite the Hamiltonian as 



Wp)|2)(2|, 



V=hA c \l)(l\ + h(A c - A p )\2)(2\ 
-fi(G|l)<3|+ fl |l)<2| +H.c.) . 



(3a) 
(3b) 



(4) 



Here, Aj for j G {c,p} are the detunings of the fields 
defined as 



A p — Lt)\2 — UJ p . 
A c = CJ13 - LO c . 



(5a) 
(5b) 



The Rabi frequencies of the probe and control fields 2g 
and 2G are related to the slowly varying amplitudes of 



£ p and £ c as 



2g = 
2G = 



2d\2 ■ £ p 

h 

2<ii3 • £c 



(6a) 
(6b) 



As shown in Fig. \V[b), the upper atomic state |1) de- 
cays to 1 3) with rate 71, and to \2) with rate 72. Further- 
more, is the decay rate of the ground state atomic 
coherence. In order to incorporate these incoherent pro- 
cesses, we use the master equation approach to describe 
the dynamics of the atomic populations and coherences 
of the three level atoms. From Eq. (jU), the equations of 
motion for the elements of the density matrix p can be 
written as 



pn = - (71 + 12) pn + WP21 

- ig*p\i - iG*p 13 , 
P22 =72/311 + ig*pi2 - WP21 , 



iGp, 



31 



Pl2 = 



Pl3 



^(71 +72) + iA p 
tGp 32 - igpn , 
7^(71 +72) + iA c 



P12 + WP22 



Pl3 + WP23 



+ iGp 33 - iGpn , 

f>23 = ~ P - l(A p - A c )]p 2 3 + W*Pl3 



iGp. 
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(7a) 
(7b) 



(7c) 



(7d) 
(7e) 



The third population follows from the conservation law 
pit + P22 + P33 = 1. In the following, we for simplicity 
restrict the analysis to the case 71 =72 = 7- 



B. Susceptibility of a thermal medium 

In the steady state, the density matrix equations ([7|1 
can be solved analytically. We are in particular inter- 
ested in the steady state value of the atomic coherence 



P12, which determines the probe field susceptibility \ of 
the medium at frequency uj p . The analytical expression 
for x to all orders in the control and probe Rabi frequen- 
cies is given in Appendix fA] The susceptibility crucially 
depends on the two photon detuning of the Raman sys- 
tem defined as 

S R = A p - A c . (8) 

For a warm atomic medium, the thermal motion of the 
atoms causes inhomogeneous broadening of the atomic 
spectra. This broadening plays a important role when 
considering the optical properties of the medium to deter- 
mine the probe beam propagation dynamics through the 
medium (38l . [3^ | . We incorporate the motion in the sus- 
ceptibility x m Eq. (|Aip by defining velocity-dependent 
detuning of a given atom 



A p (v) 



kv , 



A c {v) = A c ~kv, 



(9a) 
(9b) 



where v is the velocity of the atom. For simplicity, we 
have assumed assume k « k p « k c . The velocities of all 
atoms obey the Maxwellian distribution 



P(kv)d(kv) 



1 



_ (Ml 
e d(kv), 



v^ttD 2 

with D is the Doppler width defined by 



D = 



Mc 2 



(10) 



(11) 



at temperature T. For a warm atomic system, the aver- 
aged susceptibility is given by 



x(kv)P(kv)d(kv) . 



(12) 



C. Beam propagation dynamics 

In general, the medium susceptibility is a function of 
the position in the vapor cell due to the spatial structure 
of the applied laser fields. In order to account for the 
spatial variations, we study the propagation of the light 
beams through the medium using Maxwell's equations. 
For this, we use the slowly varying envelope approxima- 
tion, and assume that the changes along the longitudinal 
direction z of the probe envelope is very small compared 
to the changes which occur along the transverse direc- 
tions x 7 y. The beam propagation equation for the Rabi 
frequency g of the probe field under the paraxial wave 
approximation then follows as 



dg 

dz 



ic 
2lu 



A 3 ! 

dx 2 



dy 2 



g + 2nrk( X )g, (13) 



where w = ck and c is the speed of light in free space. 
The thermal motion of the atoms is included in the the 
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r (cm) 

FIG. 2: (Color online) Normalized spatial intensity profile 
of the control laser field. The solid (dashed) line shows the 
profile at the entry (exit) plane of the vapor cell, r is the 
radial distance from the beam center. The beam waist and 
the Rayleigh length of the Laguerre-Gaussian beam are chosen 
as 120 /im and 5.7 cm, respectively. 

beam propagation by considering the average response of 
the susceptibility The terms in the parentheses in 
Eq. (|13l) are responsible for the diffraction of the probe 
beam. In order to study the effect of diffraction and 
spatial dispersion on the probe beam propagation, we 
solve Eq. (|13[) numerically using a higher order split step 
operator method [401 ] . 



III. RESULTS 

A. Tailoring the medium susceptibility using a 
shaped control beam 




5 R /y 

FIG. 3: (Color online) Real part (a) and imaginary part (b) 
of the averaged susceptibility (x) in a thermal vapor obtained 
from Eq. (|12|) in the presence of the LG control beam. The 
results are shows as a function of the radial position r in the 
beam and the Raman detuning Sn/'y. The other parameters 
are chosen as G° = 17, w c — 120/im, g° — 0.2j, A p — — I7O7, 
D = 7O7, N = I0 12 atoms/cm 3 and T = 0.0017. 



In order to illustrate how the spatial structure of the 
control field enables us to guide the weak probe beam 
through warm atomic vapor, we first study the behavior 
of the susceptibility of the medium as a function of the 
spatial variation of the control field intensity. We chose 
a doughnut shaped Laguerre-Gaussian (LG) control field 
to induce a waveguide-like refractive index pattern inside 
the atomic vapor. The Rabi frequency of the LG control 
field can be written as 



G = G — — exp 



2q 



+ it 



r = \J x 2 + y 2 , 



w z — w c \ 1 + 



z- z 
zr 



q = iz R - z + z , 



= tan 



ZR 



7rw: 



A ' 



(14a) 
(14b) 

(14c) 

(14d) 
(14e) 

(14f) 



where zq and zr are the location of the beam waist and 
the Rayleigh length, respectively. 

Figure [2] shows the intensity variation of the control 
field against radial position r at the entry and exit faces 
of the medium. For this figure, we have considered a va- 
por cell of length 5 cm and filled with Rb atoms. The 
intensity of the control field approaches zero for r — > 
and maxima occur at |r| = 0.0115 cm at the entry face of 
the medium. It should be noted that there is no signifi- 
cant change in the structure of the control field envelope 
throughout the propagation through the medium. This 
is consistent with the results in the experiment by Howell 
et al |27j ] . From the point of view of an additional probe 
beam, the control field structure renders the susceptibil- 
ity inhomogeneous along the transverse direction. 

Next, we have solved Eq. (IT2|) numerically in the pres- 
ence of a continuous wave probe field to analyze the spa- 
tial inhomogeneity of the susceptibility. Figure [3] shows 
surface plots of the real and imaginary parts of (x) as a 
function of the radial position (r) and the Raman detun- 
ing (8r). For the calculation, we have fixed the single 
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FIG. 4: (Color online) Normalized intensity of the transmit- 
ted probe beam plotted against the radial distance from the 
beam center r. The solid black line shows the profile at the 
input of the vapor cell. The other two curves show the profile 
at the output after a propagation length of 5 cm. The red 
short-dashed curve is obtained without control field, whereas 
the blue long-dashed line is obtained with a LG control laser 
field. The inset shows the variation of the unnormalized out- 
put intensity of the probe beam. The initial width and Ra- 
man detuning of the off-resonant probe field are chosen as 
=48pm and 5r — —0.0157, respectively. All other param- 



eters are as in Fig[3] 



photon probe detuning to A p = — 1707, and use typical 
parameters for the 85 Rb transition: 7 = 3tc x 10 6 rad/sec, 
D = 1.33 x 10 9 rad/sec, and N = 10 12 atoms/cm 3 . It 
can be seen from Fig. 03a) that the refractive index is 
high at radial positions symmetrically located around 
r = at which the control intensity is high. Towards 
the beam center and for larger radial distances |r|, it de- 
creases. The difference between the maximum and the 
minimum of the refractive index is of order 10 -5 , which 
is in agreement with the experimental results of Howell 
et. al 27]. This shows that a LG control beam allows 
one to mimic a fibre like refractive index profile inside 
the atomic medium. From Fig. Elb) we find that the ab- 
sorption of the medium has a minimum at Sr = — O.O27, 
whereas an increase or decrease of the Raman detuning 
from this value results in increased absorption by the 
medium. 



0.75 



CD 
N 

ra 
E 
o 
z 




0.25 



r (cm) 



0.8 



0.6 



= 0.4h 

cu 

E 



0.2- 



Input beam " 


E 


i- 


! 1 1 1 1 1 ! 
— t l - 


In absence of 

control beam 1 1 

In presence of j j 

— LG control beam . . 
with w c =200|im | | / 

i i ' ' 

■ — -In presence of ! . / | 
LG control beam | |j r 
_ with w c =l00f-irn j ^ 


S f, 


■i 
i ■ 
• i 
i ■ 
■ i 

U I - 


1 - / V 'i 1 

-. .1 }L r. ; 


i! !' ii 
1 i 1 M 

j !l !,' 
! 1 !i 


1 ii • 
?! 




1 1 1 

! S 


i i 1 M 


i - 

v 




i ii 


! !!! 


11 

;> 




n i 
ii i 


i i hi 


: ' i 
i i ■ 

! 'i 


: i 


. u . 





50 c 

30. m 

20 | 
10 | 
H 



-0.05 



-0.025 





r (cm) 

(b) 



0.025 



0.05 



FIG. 5: (Color online) Normalized intensity of the transmit- 
ted probe beam. In (a), the intensity profile is shown against 
the transverse coordinate x in the y — plane at different 
propagation distances z. (b) compares the probe beam pro- 
file with and without the LG control field at the output of the 
vapor cell after a propagation of 2.5 cm. The variation of the 
unnormalized output intensity of the probe beam is shown in 
the inset. The initial amplitude and width of the off-resonant 
control and probe fields are chosen as G° = 0.757, 9° = 0.2y, 
w p — 48fim, w c — 100/im and 200/^m. All other parameters 
are the same as in Fig [3] 



1, Single Gaussian probe field mode 



B. Probe beam propagation in the tailored medium 

The analytical expressions Eqs. (fT4|) for the control 
field Rabi frequency allow us to study the effect of the 
spatial structure of the control beam on the probe beam 
propagation dynamics. For this, we have included the 
control field profile Eqs. (|14[) in the propagation equa- 
tion IP]). 



First, we study the propagation dynamics of a Gaus- 
sian probe beam. The shape of the probe beam is defined 
as 

g = g e -r*/vi(z)\ (15a) 



w 1 (z) = w p \ 1+ — , (15b) 
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FIG. 6: (Color online) Intensity of the transmitted probe 
beam in the x — y plane transversal to the beam propaga- 
tion direction at the output of the vapor cell in the absence 
(upper panel) and presence (lower panel) of the LG control 
beam, respectively. The parameters are as in Fig. [S] 



where W\(z) is the beam width. Numerical results are 
shown in Fig. 01 for a Rayleigh length of the probe beam 
chosen as 0.9 cm. The figure compares the propagation 
of the probe beam through the medium in the presence 
and the absence of the control field. In free space or in 
the absence of the control beam, we observe from Fig. 0] 
that the probe beam spreads to nearly 5.6Aw p during 
the propagation over a distance of 5 cm. The most in- 
tuitive way to understand the spreading of the beam 
is by decomposing it into a set of plane waves. In the 
course of the propagation, each plane wave acquires its 
own unique phase shift. The superposition of the phase 
shifted plane waves at the output gives rise to spread- 
ing. In contrast, for Bessel, Mathieu and Airy beams 
each constituent plane wave acquires exactly the same 
phase shift. In these cases, diffractionless propagation is 



possible through free space. 

Analogously, the manipulation of diffraction-induced 
spreading is possible by changing the spatial dispersion of 
the medium. This spatial dispersion can be substantially 
modified by employing a suitable structure of the control 
field. In particular, a suitable structure of the control 
field can cancel or even reverse the effect of diffraction. 
Our numerical results in Fig. H] obtained with the LG 
control field show that the width of the probe beam is 
37/im after a 5 cm propagation distance. Thus, the width 
of the output probe beam in the presence of the LG con- 
trol beam is narrower than the width of the input probe 
beam, demonstrating that the LG control field can be 
used to reverse and even over-compensate the effect of 
diffraction. These results are in good agreement with 
the experiment of Howell et al [I?} . We also notice from 
the inset of Fig. [4] that the integrated transmission inten- 
sity of the output probe beam is 44% of the integrated 
intensity of the input probe beam. This indicates that 
the probe beam can be guided efficiently in an atomic 
vapor. 



2. Double Gaussian probe field mode 

Next, we study the propagation of a beam with two 
Gaussian modes, using the same method as for the single 
Gaussian probe beam. Figure 5(a) shows the normal- 



ized intensity of both transmitted Gaussian beams in the 
y = plane along the z axis for a propagation length 
up to 2.5 cm in the presence of the LG control field. In 
comparison, without the control beam or in free space 
the incident probe beam undergoes sprea ding due to its 
geometrical shape as shown in Fig. |5(b)| It can be seen 
from Fig |5(b)| that again, the spreading due to diffraction 
can be eliminated using the LG control beam. The inset 
of Fig. |5(b)| shows the unnormalized transmission of the 
two Gaussian beams in the y = plane. It is clear from 
the inset of the Fig. 5(b) that the transmitted energy 
of probe beam gets reduced by 50% while the width of 
the LG control beam is changed from w c — 100 /im to 
200 /xm. We also notice that the spacing between two 
peaks of the output probe beam can be increased by de- 
creasing the width of the LG control beam. This shift in 
the probe beam peak positions arises if the probe beam 
spectrum is located on the sides of the control field in- 
tensity pattern, rather than symmetrically around the 
control field intensity maximum. As an application, the 
finesse of the output probe beam can be controlled by 
changing the width of the LG control beam. 

Figure [5] shows the radial distribution of the output 
probe beams at z = 2.5 cm, for a medium without and 



with the LG control beam, respectively. As in Fig. 5(b) 
it can be seen that in free space or without the control 
beam, the beam profiles are severely distorted. In con- 
trast, with the LG control beam, the shape of the two 
Gaussian beams is preserved. This shows that the spa- 
tial structure of the control field can be used to transmit 
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FIG. 7: (Color online) (a) Image of a non-Gaussian multi- 
peak beam at the output of a vapor cell with length L=2.5 
cm. (b) shows a comparison of the cases with and without 
control beam at the output of the vapor cell for the non- 
Gaussian beam. The parameters are as in Fig. [5] expect that 
the non-Gaussian probe beam has width 35 jUm and the LG 
control beam has width 200 /jm. 



more complicated structures over a large distance with- 
out much diffraction and attenuation. 



3. Non Gaussian probe field mode 

We proceed with a probe beam formed out of multi- 
ple secant-hyperbolically shaped patterns as an example 
for the transmission of a non-Gaussian field mode. The 
probe beam is denned as 



7° sech 



i=0 



yj{x - Xj) 2 + y 2 



(16) 



where n, Xi and Wi are the total number, the center, and 
the widths of the peaks, respectively. For simplicity, in 
our calculation, we assume that the individual peaks have 



of probe beam intensity in the presence of the LG con- 
trol beam at the exit plane of the vapor at z — 2.5 cm. 
A corresponding comparison of the non-Gaussian beam 
propagation through the medium with and without the 
control beam is shown in Fig. |7(b)| It is important to 
note from Fig. |7(b)| that the output beam retains sub- 
peaks of secant-hyperbolic profile as they were present 
in the input beam. Thus the LG control beam allows 
to propagate Gaussian as well as non-Gausssian beams 
through the medium without diffraction. 



IV. SUMMARY 

We have discussed diffraction-less propagation of 
Gaussian and non-Gaussian probe beams through a 
warm atomic vapor suitably tailored using an additional 
control beam. The Laguerre-Gaussian control beam 
writes a spatially modulated index of refraction in the 
vapor, which leads to a modified dispersion experienced 
by the probe beam, and subsequently to diffraction- 
less propagation of the probe beam. We have shown 
that probe beams of single Gaussian, multiple Gaussion, 
and multiple non-Gaussian structure can be propagated 
through the medium without diffraction. 
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Appendix A: Steady state susceptibility 

The steady-state susceptibility is given by 
N|d| 2 Af 



X 



h V 



(Al) 



where the numerator Af and the denominator T> evaluate 
to 

Af =\G\ 2 [ 7 (*T - 8 R ) { |G| 2 + ( 7 - *A p )(r - i5 R )} 

+ \g\ 2 { 7 (*r - 6 R ) + T(A C + A p )}] , (A2) 

D -7lff| 6 + M 4 {3|G| 2 ( 7 + 2r) + 2 7 ( 7 r + 6 R A C )} 

+ 7 |G| 2 {(r 2 + 4)(7 2 + A 2 ) + 2|Gf ( 7 r + S R A P ) 
+ |G| 4 } + | 5 | 2 {3|G| 4 ( 7 + 2r) 

+ 7 ( 7 2 + A 2 ) (r 2 + 5\) + (4 7 2 r + 7 (er 2 + 44) 

+r(A c + A p ) 2 )|G| 2 } . (A3) 
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